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ABSTRACT 

It is interesting to superimpose the Pauli-Villars regularization on the lat- 
tice regularization. We illustrate how this scheme works by evaluating the axial 
anomaly in a simple lattice fermion model, the Pauli-Villars Lagrangian with a 
gauge non-invariant Wilson term. The gauge non-invariance of the axial anomaly, 
caused by the Wilson term, is remedied by a compensation among Pauli-Villars 
regulators in the continuum limit. A subtlety in Frolov-Slavnov's scheme for an 
odd number of chiral fermions in an anomaly free complex gauge representation, 
which requires an infinite number of regulators, is briefly mentioned. 
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It seems interesting to put a Pauli-Villars type Lagrangian level regulariza- 
tion on the lattice. The interest is twofold: The Pauli-Villars regularization [1] 
for fermion one-loop diagrams can be expressed as a Lagrangian of regulators 
(bosonic and fermionic spinors). In actual perturbative calculations however, the 
Lagrangian has to be supplemented with additional prescriptions, such that the 
momentum of propagators have to be assigned in the same way for all the fields, 
and the integrand in the momentum integral has to be summed before the inte- 
gration. Once the Lagrangian is put on the lattice, no prescription are needed and 
one is free to choose any momentum assignment. (To get a finite gauge invariant 
result in the continuum limit one has to assign the momenta of all fields on the 
lattice in the same way.) 

More interestingly and importantly, "superimposing" a different kind of regu- 
larization on the lattice regularization may give some clue to the lattice regulariza- 
tion of chiral gauge theories. No manifestly gauge invariant lattice formulation of 
the chiral gauge theory, being consistent with the unitarity and the locality, is yet 
known [2]. In particular, for a chiral fermion in a complex gauge representation, it 
is impossible to introduce in a gauge invariant way the Wilson term [3] to eliminate 
unwanted species doublers.^ The difficulty of a manifestly gauge invariant lattice 
formulation of chiral gauge theories is highlighted by the No-Go theorem [4]. 

The basic idea of "superimposing" is quite simple. Let us consider, for example, 
the naive momentum cutoff regularization applied to fermion one-loop diagrams 
in QED. This regularization breaks the gauge invariance, generating gauge non- 
invariant contributions. However we may use in addition say, the gauge invariant 
dimensional regularization. With this superimposed regularization, the infinite 
momentum cutoff limit can be taken and we are left with gauge invariant expres- 
sions in the dimensional regularization. Of course there is no real need to break the 
gauge invariance by introducing the momentum cutoff in this example, but with the 

f For chiral fermions in a real-positive gauge representation, and for even number of chi- 
ral fermions in a pscudoreal representation, it is possible to introduce a gauge invariant 
(Majorana-type) Wilson term. 
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lattice regularization, it is not obvious how to treat chiral fermions in a manifestly 
gauge invariant manner. To perform this program it is clearly crucial that there 
exists a regularization which preserves the gauge symmetry and simultaneously is 
congenial to the lattice regularization. 

In fact a proposal based on this idea has been made by Frolov and Slavnov [5] 
(see also [6]). They used the gauge invariant generalized Pauli-Villars regular- 
ization [7,8] for chiral fermions in an anomaly free complex representation, and 
discussed that taking the continuum limit a — > (a is the lattice spacing) with 
an appropriately scaled regulator mass M(a) Cl/a, the regulator fields compen- 
sate the effect of gauge non-invariant Wilson term and that the gauge invariant 
regularized continuum theory [7] is reproduced. 

We illustrate in this letter how this scheme works by evaluating the axial 
anomaly [9,10,11] in a simpler lattice fermion model, the Pauli-Villars Lagrangian 
with a gauge non-invariant Wilson term in the lattice vector gauge theory. The 
relation to the proposal in [5] and the possible implication will be commented on 
later. 

Before considering the Pauli-Villars regulators, let us study for a while a mas- 
sive Dirac fermion coupled to a background gauge field and the axial U(l) Ward 
identity. The naive fermion action is 

1 = ^2^(x)[i]p(x) -m]i/>(x) = ^2^(x)[-ilp{x) -m]ifj(x), (1) 

X X 

where the covariant derivative on the lattice D^(x) has been defined by 
D^x) = J- \u»{x)e ad » - e- a ^uUx) 



2a 

D,(x) ee -i- \u,(x)e- a ^ - e^U^x) 



(2) 



In the continuum limit a — * 0, we parameterize the link variable as U^x) = 
e ta 9 A ^( x ). As is well known, the naive action (1) contains unwanted species dou- 
blers. Therefore we add the Wilson term to decouple them in the continuum limit, 
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but an artificially chosen gauge non-invariant one: 



Iw = ^(x)R(x)^(x) = ^VW(-SW)^), (3) 



with 



r w = ^E( e " a " + e ~ ad * - 2 )' r( x ) = -^J2{ e ' a ^ + ^ - 2 )- ( 4 ) 

Although iyj/ is irrelevant in the naive continuum limit, the effect of hard breaking 
of chiral and gauge symmetries survives in the axial anomaly as we will see below. 

The axial U(l) Ward identity for the lattice action / + I\y is derived by per- 
forming a change of variable ip(x) — > e ta ( x ^ 5 if}(x) and ip(x) — > ip(x)e ia ^ x ^ 5 in the 
partition function. For an infinitesimal a(x), the action changes as 



X 

I w -> Iw -^2a(x)B(x), 

X 

where the divergence of the axial U(l) current has been defined by 

d^(x) = tp(x)lp(x)~f 5 ip(x) + tp(x) $(x)^(x) 
= Y,d i ^{x) 1 ^{x)\ +0(a), 



(5) 



(6) 



(the second line is the naive continuum limit) and 

B(x) = -ijj(x)i^ 5 R(x)ijj(x) + ip{x)i'y 5 R{x)ip(x), (7) 

the explicit axial U(l) breaking part. Since the functional integration measure is 
invariant under the change of variable with the lattice regularization, the Ward 
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identity reads 

(Vs 0*0) = (^(x) l5 ^(x)) + (B(x)) . (8) 

Let us evaluate the right hand side of (8). We first concentrate on the vacuum 
expectation value of B(x)* 



(B(x)) = tri^R(x) (ijj(x)i)(y))\ x=y -tv(ijj(x)ij(y))i-f 5 R(y)\ x=y 
= b(x) +*b(x). 



(9) 



They are evaluated by the lattice propagator in the presence of background gauge 
field, 

b(x) = — tr i'jKR(x) -—8(x,y)\ , 

K ' lb V ! 0{x) -m + R(x) V 

1 - (10) 

b{x) = -trS(x,y)— — il5 R(y)\ x=y , 

i0(y)+m + R(y) 

where the delta function on the lattice is defined by 5(x, y) = 5 x ^ y / a 4 = J^^ a d 4 k x 
e i*(a:-y)/(27r) 4 and hence 



A 



b(x) 



x 



^e- tkx / a t l5 R(-zlp-m + R) 



-E,Dl + (m- Rf + £^[ 7 M 7"] [Dp, A,]/4 + R] 



ikx/a 



In deriving the above expression, we have multiplied (—iJp — m + R) on the numer- 
ator and on the denominator, and used a relation Jj) 2 = — ^2 D 2 + ^ ^[7^, 7"] x 
[Dp, A,]/4. Next noting 



e -ikx/a D ikx/a = 1 gin k + n e -ikx/a Re ikx/a = r V- (cog k X ) + Rj ( 12 ) 



★ Our calculation method is similar to that of [11], but seems rather simpler. 
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where 



and 



D, 



2a 



e ik *(U^e ad » - 1) + e _i *"(l - e'^lA) 



cos kn(d^ + igAn) + 0(a), 



R = 



2a ^ L 



l)-e-^(l-e- a ^) 



J^sin fc^ + 0(a), 



(13) 



(14) 



we have 



b(x) 



j ,f [y y^(c// - 1) + a-R ^ + r y^(c^ - 1) - - am + aR 



x |- y^(jg p + a_D p ) 2 + ^y^(c p - 1) - am + ai? 



p,a 



(15) 

where the trigonometric functions have been abbreviated as = sin fc M and c M = 

cos 

In (15), the expansion with respect to a is straightforward because the trace of 
gamma matrices requires at least four of them (tr 757^7^7^7°" = — Ae^ vpa ). Finally 
using 

[Dp, D v \ = igc^CyF^y + 0(a), 

[D^, R]=rgJ2 W„{d v Ap) + 0(a), (16) 



where F^ = d^A u - d v A il + ig[A^ A v \, 



Jim b(x) = -j^h^e^trF^ - yJ^W^tr dfl A„F p(7 . (17) 
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In (17), h(r) and hij') are the well-known lattice integrals [9,11]: 




h(r) 




[Ed -<*)]'} 

1 p p ) 



(18) 



Hr) 



P v=£p 



-IB 1 -<*>]'} 

k p p } 



3 ' 



and obey h(r) + h(r) = — 7T 2 /2. A similar calculation shows b(x) = b(x). Thus 
finally we arrive at 



Several comments are in order: It is well-known that, with the gauge invariant 
Wilson term [3], the species doublers effectively act as the Pauli-Villars regulators 
(i.e., with an alternative axial charge [9]) and the gauge invariant correct axial 
anomaly is reproduced in the continuum limit [9,11]. In our present case, the last 
term on the right hand side of (19) is not gauge invariant (the first term will be 
shown to be gauge invariant). It would be gauge invariant if the last coefficient 
was 2/3 instead of 1. It also depends on the Wilson parameter r, although for an 
infinitesimal r the term vanishes lim r ^o-^2( r ) — 0. Therefore the effect of hard 
breaking of chiral and gauge symmetries in the Wilson term (3) survives in the 
continuum limit. In other words, our identification of the axial current (6) was not 
gauge invariant (with the gauge non- invariant Wilson term (3)), and the operator 
does not coincide with the naive continuum limit. We may redefine the axial current 
as only the first line in (19) survives and to restore the gauge invariance. Such an 



lim (d„J^(x 



)> 




(19) 
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intricacy of an operator identification in the lattice regularization in the view point 
of the anomaly and the usefulness of superimposing a different regularization are 
emphasized in [12]. 

For our purpose, the first term in (19) is important: 



(2imip(x)^5ip(x)) 



= 2i tr —r 



d 4 k 



a 4 J (2tt 



£ + r ^^(c^ — 1) — ialp — am + ai?j 



x 



|— ^^(is p + aDp) 2 + \r '5~](cp — 1) — am + aR 



(20) 

A simple expansion by a is impossible in (20) because of the singular infrared 
behavior near k ~ 0. Thus we divide the integration region to the "outer" region 
\k\ > 5, which is free of the infrared divergence, and the "inner" region \k\ < 5 
with an infinitesimal (5 C 1 [9], In the outer region, we may safely expand the 
integrand with respect to a, yielding 



lim {2im^{x) 1 ^{x)) mtex = 0. 



(21) 



In the inner region we expand it by the external gauge field: 



<2im^(x)75^(a;)) iriner = 2imT 5 

oo 1 n r „ ,4 

n=\ i=l *-Xi,Hi J \ J 



iapi/2 



rr 2 '" /1 "(PbP2,--',Pn). 

(22) 
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For example, Ts = Tg (p) = 0, and 



rf(p,g) 



= - tr 75 



5/a 



d 4 k r 



-5/a 



(2tt)< 



S(k+p + q)V»{k +p + q,k + q)S(k + q)V l '{k + q, k)S(k) 
+ (pL <-> v,p q) 



(23) 



where 



S{k) = s i n + m H — ^^(1 — cos a/c M ) , 



7 M cos a 



-/ci„ H — ko,?] + f sin a( -ki,, H — fco,,^ 



(24) 



Since what to be worried about is the infrared divergence, we may simply expand 
the numerator by a as V^(k + p + q, k + q) = + 0(a) etc. In the denominator, 
we may expand the propagator as, S(k+p) = E^^+Pm) + m + 0(5 2 )}~ 1 because 
5 <C 1. Note that the subleading contributions from the denominator always give 
a ultraviolet convergent integral. From these arguments, we have 



lin l r r(P'?) 



5/a 



—g tr75 lim 



d 4 k 



x 



a"b J (27r) 4 

-5/a 

1 



4tt 2 



k'+p'+ti + m ^ + i + m # + 
l l 

m^e^ po p p q a J 2ydy J dx 



m 



(fi <-> <-> g) 



+ 0(5 2 ) (25) 



X 



+ 0(5 2 ). 



m 2 — y(l — y)g 2 — 2xy(l — y)p • q — xy{l — xy)p 2 
After the safe limit 5 — > 0, (25) is nothing but the expression in the continuum 



★ The so-called anomalous vertex [9] docs not contribute to this function. 
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theory. The same consideration can be repeated for higher point functions, and we 
may summarize this fact compactly as 



lim (2im^(x)^5ip(x)) = 2i lim trm75— r 5(x — y), (26) 

a^o y^x iJp c — m 

where Tj) c = 7 M (<9^ + igA^) is the covariant derivative in the continuum theory. 
Note (26) is finite without any further regularization as the last line of (25) shows. 
The first term in (19) is therefore gauge invariant. As we have seen in (21), only the 
physical fermion near k ~ contributes to the operator (20), and the effect of Wil- 
son term which couples only to the species doublers, is invisible in the continuum 
limit. 

Combining (19) and (26), we have for a single massive Dirac fermion, 
= 2i^tTmTs m -—6(x-y) (27) 



y->x i]p c — m 

• 2 -2 

+ ^n^^F^ + ^I 2 (r)e^tr[d^A 1/ d p A (7 + igA v A p A (T )}. 



Let us now introduce the Pauli-Villars Lagrangian on the lattice. We introduce 
the regulator fields ip n , where n — 1, 2, • • • , N, and i/jq the original fermion to be 
regularized. We assign the even number index for fermionic fields and the odd 
index for bosonic ones. We also denote the mass of those fields as m n and assume 
the masses of the regulator fields are of the order of the "cutoff" parameter A. The 
Pauli-Villars regularization condition [1] requires J2n=o(~ 1)™ = J2n=o(~ l) m n = 
0. For the combined system, we can use (27) for each fermionic field and that with 
a reversed sign for each bosonic field. For example, by summing up the first term, 
we have 

N 



2i lim tr V(-l) n m n75 — S(x - y) 

= <22m ^o75^0> + 2t lim tr l5 f(lp 2 c /A 2 )5(x - y), 

v ' y— >x 



(28) 
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where we have defined the regulator function 

m . f (29) 

It follows from the definition and the Pauli-Villars condition, f(0) = 1 and f(t) = 
rriQ/(A 2 t) + 0(l/t 2 ). For example, a possible choice is m\ = A 2 , m\ = 2A 2 , 
m§ = A 2 + m\, and lim A -oo f(t) = 2/(t + l)(t + 2). 

On the other hand, the second line in (27) cancels out among the fermionic 
and the bosonic fields because it is independent of the mass of the field. Therefore 
the total axial anomaly is given by 

lim hm ld^ix)) = (2tm<rt 0l ^ ) + ^e^ pa tr F^F pa . (30) 

A— >oo a^O 107T 

In deriving this, we have used the fact /(0) = 1 and /(oo) = f'(oo) = • ■ ■ = 
for A — > oo and the actual calculation is identical to that of [13]. We have thus 
recovered the correct gauge invariant form of the axial anomaly. The present 
analysis may be repeated for the conformal anomaly of the Wilson fermion, for 
which the Wilson term gives rise —15 times the correct coefficient [12]. We expect 
that the correct coefficient will be reproduced in the present scheme, because the 
real Pauli-Villars regulators, i.e., with an alternative statistics, always eliminate 
the effect of Wilson term. 

Finally we briefly comment on the implication of above demonstration for the 
proposal in [5] . By a suitable change of variable, 

1 



X (x) = -j= [Pr^{x) + P L CT n C D ^ T (x) 
<p(x) = ±= [Pr<P(x) + P L CC D f{x) 



(31) 



P R ct>{x) + P L CC D ^{x) 



(x{ x ) is a fermionic field and <p(x), <p(x) are bosonic fields), it is possible to rewrite 
the lattice Pauli-Villars Lagrangian of [5] basically in the form of / + Iw in (1) 
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and (3). Note however that the number of the degree of freedom is doubled in (31) 
because, among four chirality components of x( x ) an d X.{ x )i only two of them are 
independent (Pr?P(x) and tp(x)Pi are independent in the original chiral model). 
Therefore the above result (27) divided by two with m = might be regarded as 
the fermion number U(l) anomaly [14] of the original massless chiral fermion in 
the lattice model of [5] . 

/,From the analyses in [7,8], we know that it is possible to construct a gauge 
invariant Pauli-Villars Lagrangian, by utilizing a finite number of regulator fields, 
for chiral fermions in a real-positive representation, and for an even number of 
chiral fermions in a pseudoreal or an anomaly free complex representation. In 
these cases the situation would be the same as the above analysis and we would 
have (30) with mo = as the twice of the fermion number anomaly, i.e., the correct 
result, lim A ^oo lim a ^o (Vfa)) = W 2 e^ pa tr F IIV F pa /(32n 2 ). 

For an odd number of chiral fermions in a pseudoreal representation and in 
an anomaly free complex representation, it is necessary to introduce an infinite 
number of regulators [7,8]. For the former case, it is possible to first put a finite 
number of them on the lattice, in a way that Yln=o(~ -0™ = 0> anc ^ ^hen to ta ^ e 
the limit N — > oo. Thus the above result might hold even in the N — > oo limit. 
On the other hand, the situation seems more subtle for the latter case, because 
the condition ^2n=o(~^) n = is never satisfied for a finite iV [7,5]. The reason 
is that, while the original chiral fermion belongs to a complex representation, all 
the regulator fields belong to a "doubled" representation [7,5,8], namely the con- 
tribution of one regulator field is twice of the original fermion. For example, if 
we put the same finite number of fermionic and bosonic regulators first, and then 
take the limit N — + oo, it would correspond to lini/v->oo ^2n=o(~ -0 n = 1 i n the 
above notation. On the other hand, the regulator function (29) would be given by 
limjv-^oo f(t) = n\fi/ sinh(7r\/t) [5]. Therefore the fermion number anomaly would 
be given by 
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(32) 



and the gauge non-invariant piece survives. The conclusion in [5], on the other 
hand, would imply 



The above (admittedly handwaving) argument, (32) and (33), shows that the 
anomaly is quite sensitive on detailed way of the limit a — > and iV — » oo, in 
Frolov-Slavnov's scheme for an odd number of chiral fermions in an anomaly free 
complex representation. 

In conclusion, we have illustrated how the superimposing of the Pauli-Villars 
regularization on the lattice regularization works, utilizing the axial U(l) identity. 
Simultaneously, we expect that the scheme also improves non-anomalous chiral 
symmetric properties of the Wilson fermion in QCD. To pursue this program fur- 
ther, however, we have to treat for example, the fermion self-energy, for which the 
Pauli-Villars regularization gives no clue. On this point, a superimposing of the 
higher covariant derivative regularization has been proposed [5]. 

We thank S. Kanno for helpful discussions. The work of H.S. is supported in 
part by the Ministry of Education, Science, Sports and Culture Grant-in-Aid for 
Scientific Research, No. 08240207, No. 08640347 and No. 07304029. 




im (V(a;)) 




(33) 
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